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Abstract — A measure for the maximum quantum information 
transfer capacity (ITC) between nodes of a spin network is 
defined, and sliown to induce a metric on a space of equivalence 
classes of nodes for liomogeneous chains with XX and Heisenberg 
couplings. The geometry and curvature of spin chains with 
respect of this metric are studied and compared to the physical 
network geometry. For general networks hierarchical clustering 
is used to elucidate the proximity of nodes with regard to the 
maximum ITC. Finally, it is shown how minimal control can be 
used to overcome intrinsic limitations and speed up information 
transfer. 

L Introduction 

Networks of interacting quantum particles — so-called spin- 
networks — are important for transferring and distributing 
quantum information between different parts of a larger system 
such as different quantum components on a chip ||T]. Spin 
chains, linear arrangements of spins, for example, can play 
the role of classical wkes connecting two parts, and branched 
networks allow the distribution of quantum information to 
different nodes. The way quantum information propagates 
through spin networks, however, is quite different from clas- 
sical information flow due to quantum inference effects. In 
particular, quantum state transfer between the nodes of the net- 
work is limited by fundamental physical principles and perfect 
quantum state transfer is usually possible only in very special 
cases. Quantifying state transfer fidelities for spin networks 
is not easy in general, but for certain types of networks such 
as spin-^ particles with interactions of so-called XXZ type, 
for example, this problem can be reduced to the maximum 
probability for a single excitation to propagate from one node 
to another, which can be computed efficiently numerically, 
and in some cases analytically. Thus, for the respective spin 
networks, the latter is a basic measure for the maximum 
information transfer capacity (ITC) between different nodes. 
For certain types of simple networks such as homogeneous 
chains the maximum ITC is shown to induce a metric on a 
set of equivalence classes. We study the topology of simple 
networks with regard to this metric, showing that it differs 
substantially from the physical geometry. Networks with trivial 
physical geometry such as linear chains can have a surprisingly 
rich geometric structure including curvature. Analysis of the 
latter shows that spin chains appear to be Gromov-hyperbolic 
with regard to this maximum ITC metric but unlike classical 
hyperbolic networks their Gromov boundary appears to be a 



single point. For more complex networks the maximum ITC 
does not induce a metric but we can use hierarchical clustering 
to assess the proximity of nodes with regard to information 
transfer capacity. Again, the resulting cluster structures differ 
substantially from the neighborhood relations induced by the 
physical geometry of the network, showing that the latter is not 
very useful in assessing the maximum ITC between nodes in a 
network, unlike in the classical case. Finally, we consider how 
minimal control of a single node in the network allows us to 
change information flow in the network, effectively changing 
the information transfer capacities (and network "topology"), 
allowing us to achieve higher state transfer fidelities as well 
as generally speeding up the rate of information transfer. 

II. Information Transfer Capacity 

The Hilbert space of spin networks with XXZ interactions 
can be decomposed into so-called excitation subspaces H = 
®n=o^n, where n is the number of excitations in the network, 
ranging from to N . If we denote the spin basis states by |t) 
and 14^), taking the latter to denote the ground state, then the 
0-excitation subspace consists of a single state |0) := \\r) ® 
■ ■ ■ (g) 14,), while the one-excitation subspace consists of 
states \n) :— \\) ® ■ ■ ■ ® \\) ® ■ ■ ■ ® \\), where the excitation 
It) is in the nth position. Thus transferring a quantum state 
I'i/') = cos(6')||) +e"^ sin(6')|t) from spin m to n is equivalent 
to transferring an excitation |t) from spin m to n: 

|^)„ = cos(0)|O) +e'^ sin(0)|m) 

^ cos(6')|0) + e*-^ sin(e')ln) |i/;)„, 

where we used the shorthand |i/')„ to denote a product state 
\\r) ® I^A) ® li), whose nth factor is all others being jl). 

The probability that an excitation created at site \i) has 
propagated to site |j) after some time t is given by 

p(N>,|j),i)-|01e-'^*N)P (1) 

in the system of units where h ~ 1. The maximum of 
this probability p{i,j) = maxt>op{i,j,t), or a monotonic 
function thereof such as 

d{i,j) = -'iogp{i,j) (2) 

gives a measure for the maximum state transfer fidelity be- 
tween two nodes in a spin network without control, quantifying 
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TABLE I: ITC geometry of linear Heisenberg chains. 




Chain length N Chain length N 



Fig. 1: Chain length versus 4-point Gromov 6 (left) and scaled 4-point Gromov S (right). Scaled-Gromov S remains below 
the upper bound singled out in |4J, 4-point Gromov S saturates at large chain length, revealing Gromov hyperbolic property. 
G ^ L + M + S where L, M and S are the pairs of opposite diagonals of a quadrangle corresponding to the largest, medium 
and smallest length. 



the intrinsic capacity of a spin network for quantum state 
transfer tasks. It can be shown that 



we 
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where H — ^fe kfe) {"^k I is the eigendecomposition of the 

Hamiltonian operator H of the network in the first excitation 
subspace Hi. If the rescaled eigenvalues are 
rationally independent |j2] Proposition 1.4.1] then it can be 
shown that for any e > there exists a > such that 
Pmaxihj) — p{h j,t<i)) < i-C-i the maximum ITC is a tight 
bound and attainable in the limit. 

III. ITC Geometry & Curvature of Chains 

Inspired by |3| we may hope that the maximum ITC 
measure d{i,j) defined above can be shown to be a distance. 
d{i,j) clearly satisfies d{i,i) = as p{i,i) = 1, and we 
also have symmetry as obviously p{i,j) = p{j,i)- However, 
d{i,j) can vanish for i 7^ j if p{i,j) — 1, and in general we 
cannot expect the triangle inequality to hold. In special cases. 



however, such as for homogeneous chains with either XX or 
Heisenberg coupling, numerical exploration for systems up to 
500 spins shows that the triangle inequality 



Ajfe = k) + d{j, k)) - d{i,j) > 0. 



(3) 



seems to be universally satisfied. This renders d{i,j) a semi- 
distance, which induces a proper distance on a set of equiv- 
alence classes defined by identifying nodes with d{i,j) — 0. 
Specifically, for Heisenberg or XX chain with uniform cou- 
pling the N spins can be shown to form [^] equivalence 
classes comprised of spins n and TV + 1 — n, which we shall 
denote by n for n = 1, . . . , f-^^/S] in a slight abuse of notation. 
On this set of equivalence classes the ITC measure d{i,j) is 
a metric, and it is interesting to study the induced geometry 
of a spin chain with respect to the ITC metric, and how it 
differs from the physical network geometry. Table |l] gives the 
ITC geometry for Heisenberg chains up to TV = 10, showing 
that it is very different from the (trivial) physical geometry. 
A Heisenberg chain of length N ~ 7, for instance, appears 
as a pyramid structure with an equilateral triangle as its base 
formed by the vertices a — {1, 7}, b — {2, 6}, c = {3, 5} and 
d = {4} as its apex, where n refers to the index of the spin 
in the chain. Chains up to length = 8 can be embedded 




Spin index n<N Spin index n<N 

Fig. 2: Inertia (a = 2) versus chain length plots for Heisenberg chains. For N odd (left) the central node is a strong anti-gravity 
center. For N even (right) there are weaker anti-gravity centers between end spins and middle (right). The inertia is symmetric 
about the central node as antipodal spins belong to the same equivalence class. 




Spin index n<N Spin index n<N 

Fig. 3: Inertia (a = 2) versus chain length plots for XX chains. For N odd (left) the central node is a strong anti-gravity center 
and weaker anti-gravity centers exist for intermediate nodes. For N (even) the anti-gravity centers are located between end and 
central spins. Again, the inertia is symmetric about the central node as antipodal spins belong to the same equivalence class. 



in M.^ but for longer chains higher-dimensional spaces are 
required. By Schoenberg's theorem |j5J the distances d{i,j) 
can be realized in M.'^ if and only if the Gram matrix G = (Gy ) 
is positive semi-definite of rank d, where 

G,, = i {d{i, Nf + d{j, Nf - d{t, jf) . (4) 

(This is equivalent to the Cayley-Menger matrix criterion 
of Q.) We numerically verified that the Gram matrix for both 
Heisenberg and XX chains is positive semi-definite for spins 
up to = 500. Analysis of the rank of the Gram matrix 
furthermore suggests that the dimension d required to embed 
a chain of length N is ['^^j^l for both Heisenberg and XX 
chains. Note that Schoenberg's theorem also applies to the 
geometric realisability of general spin networks and whether 



our ITC measure fulfills the triangle inequality. 

To better understand the geometry for very long chains we 
analyze its curvature a la Gromov. From both the Gromov 
and the scaled-Gromov point of view ||7j spin chain of both 
Heisenberg and XX coupling type appear hyperbolic. The 
strict Gromov point of view is displayed in Fig. JTJa) and 
the scaled-Gromov point of view is shown in Fig. fUb). The 
Gromov property might appear to conflict with the Euclidean 
embeddability of the metric space made up by the clusters 
with the ITC metric, as the Bonk-Schramm |8| theorem says 
that Gromov negatively curved spaces are embeddable in hy- 
perbolic space. However, some metric spaces are embeddable 
in both Euclidean and hyperbolic spaces, the most striking 
example being that of a complete graph with uniform edge 



weight, and spin chains with uniform coupUngs appear to fall 
into this category. 

Classical communication networks, both wired j?), Q and 
wireless |[3), have been shown to be Gromov hyperbohc. 
Gromov hyperbolic spaces have a unique vertex that achieves 
the minimum inertia |10|, the so-called gravity center 

5 = argmin/(i) = argmin > c?"(i,j), a > 1. (5) 

j 

Classical networks indeed show a point of minimum inertia, 
which can be interpreted as a congestion point {\T\. Classical 
networks also have the property that their Gromov boundary 
in the asymptotic limit is a Cantor set p2) . Quite surprisingly, 
quantum networks differ from their classical counterparts in 
that they have points of maximum inertia, or anti-gravity 
centers, as shown in Figs |2] and [3] for Heisenberg and XX 
chains, respectively, and their Gromov boundary is a single 
point. The interpretation of the antigravity center is that the 
information flow in the network avoids this node, making it 
difficult to transfer excitation to and from it. It is interesting 
to note in this context that the distance between antipodal 
nodes is 0, meaning that we can achieve state transfer fidelities 
arbitrarily close to 1, i.e., near perfect state transfer between 
the ends of the chain of any length given sufficient time, while 
near perfect state transfer is never possible between any pair 
of nodes with d{i,j) > 0, no matter how long we wait. 

IV. ITC "Topology" Of General Networks 

The maximum transfer probability is also a useful measure 
for the maximum fidelity of quantum state transfer in general 
spin networks. Although it is in general no longer a metric, it 
can still be used as a similarity measure to define a hierarchical 
clustering. We use the pairwise clustering algorithm introduced 
in p3) . Pairs of nodes are grouped hierarchically into clusters 
in order of their similarity, and only those clusters whose 
elements are closer to each other than any element outside 
the cluster are preserved. If we define a relation 

\a)^,\b):^d{\a),\b))<e (6) 

for nodes |a), \b), then these clusters are the equivalence 
classes of =e for a certain e. The cluster hierarchy reveals 
the closeness of the nodes in terms of information transfer. 

For example, consider a network of 10 spins distributed 
in a square forming a general spin network as shown in 
Fig. |4] The positions of the spins are indicated by the blue 
dots. Taking the coupling strength J{i,j) between spins i 
and j to be inversely proportional to the cube of the physical 
distance between the nodes, we compute the Hamiltonian of 
the network, assuming XX-coupling. We then diagonalize this 
Hamiltonian and compute the maximum transfer probabilities 
p{i,j) and the associated d{i,j), which are used as input for 
our clustering algorithm. The resulting hierarchical clustering 
structure is shown in Fig.|4] Different colours indicate clusters 
for different similarity levels; i.e. for clusters marked in 
the same colour there exists an e for which these are the 
equivalence classes of Again, the example shows that the 




Fig. 4: Clustering induced by maximum transfer probability 
measure for a general spin network shows that clustering 
differs from what would be expected if we clustered according 
to physical distance. 



physical distance of the spins is not a good measure of their 
proximity in an quantum information transfer fidelity sense. 
For instance, spin 9 is physically closer to 3 than any of the 
other nodes, yet the clustering indices that 9 and 3 are several 
levels removed with regard to the maximum ITC. 

V. Control of Information Transfer 

To overcome intrinsic limitations on quantum state transfer 
or speed up transfer, one can either try to engineer spin chains 
or networks with non-uniform couplings or introduce dynamic 
control to change the network topology. The idea of engi- 
neered couplings was originally proposed to achieve perfect 
state transfer between the end spin in spin chain quantum 
wires | [T5| . The analysis above shows that engineering the 
couplings is not strictly necessary. As the distance between 
the end spins with regard to the ITC metric defined above is 
zero for uniform XX and Heisenberg chains, we can achieve 
arbitrary high fidelities state transfer between the end spins 
if we wait long enough. Engineering the couplings, however, 
can speed up certain state transfer tasks such as state transfer 
between the end spins at the expense of others. 

A more flexible alternative to fixed engineered couplings 
is to apply control. For instance, suppose we would like to 
transfer an excitation from node 1 to 4 for an XX chain 
of length N = 7. Node 4 being the anti-gravity center, the 
maximum transfer probability without control is low regardless 
how long we are prepared to wait. If we are able to change 
the Hamiltonian of the network by applying some control 
perturbtation so that H — Hq + u{t)Hi we can change the 
situation even if the control is restricted to u{t) = 0, 1, and 
Hi is a local perturbation, e.g., of a single spin induced by a 
magnetic field, e.g.. Hi = ctz- Switching the control on/off 
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Fig. 5: Population of |4) for a uniform chain XX chain of 
length N — 7 under free and controlled evolution. Under free 
evolution (dashed) the population cannot exceed 0.4268 (dash- 
dot line) but control can overcome this restriction resulting in 
near perfect excitation transfer 



at times t„ induces the evolution 

Uu{tn, 0) = Un-i mod 2{tn, tn-l) ' ' ' Ui{t2, ti)Uo{ti,to) 

(7) 

where Uoitk,tk^i) = e-^'"-'"-^^"" and Ui{tk,tk-i) = 
^-i(tk-tk-i)(Ho+Hi) ^ gy optimizing the control sequence, i.e., 
in this restricted case the switching times {tk}, we can change 
the dynamics to achieve near perfect transfer of an excitation 
or quantum state to a desired target node. In 1 16] it was shown 
that applying a simple bang-bang control sequence to a single 
spin can significantly speed up quantum state transfer between 
the ends of a chain, but control also allows us to overcome 
fundamental limits imposed the maximum ITC, enabling us 
to achieve near perfect excitation transfer to the anti-gravity 
center in very short time, as shown in Fig. |5] for example. 
We can think of the control sequence as implementing an 
effective Hamiltonian Hctf defined by e"*'^""*" = f7„(t„,0) 
at time i„. This effective Hamiltonian differs from the system 
Hamiltonian Hq, as do the transition probabilities. For com- 
parison, we diagonalize the effective Hamiltonian, compute 
the associated maximum transition probabilities Pcffihj), and 
use hierarchical clustering to elucitate the proximity relations 
between spins under the original and effective Hamiltonian. 
Fig. |6] shows that the results for a particular control example. 

VI. Conclusions 

Following the recent trend of geometrization of classical 
communication networks, we have here developed the geom- 
etry of spin networks using an Information Transfer Capacity 
metric. Classical and quantum networks bear the similarity 
that they are both Gromov hyperbolic, with the difference 
that classical networks have a Cantor Gromov boundary while 
spin chains have their Gromov boundary reduced to a point. 
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Fig. 6: Clustering induced by maximum transfer probability 
for a uniform chain XX chain of length N ~ 7 without control 

(a) and clustering induced by controlled transition probability 

(b) for a bang-bang control sequence designed to achieve 
perfect state transfer from spin 1 to 4. 



This probably accounts for yet another discrepancy: classical 
networks have a gravity center (a congestion point) while 
quantum networks have an anti-gravity center, a spin diffi- 
cult to communicate with. The broader implication of this 
geometrization is that it specifies where control is necessary 
to overcome such limitations of the physics. 
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Appendix 

A. proof of attainability 

Theorem 1: If the numbers are rationally inde- 

pendent ||2j Proposition 1.4.1], then Ve > there exists a i > 
large enough such that p{i,j) — p{i,j,t) < e. 

Proof: Set the system of units such that h = 1. In order 
to reach the maximum probability exactly, one must find a t 
such that 



-Xkt ^ (2m + l)7r if Sgii{vikVjk) 



-1 



-Xkt = (2m)7r if Sgn{vikVjk) = +1 

where Vik = {i\vk) is the projection of the eigenstate Vk onto 
the basis state etc, and — Vik as the eigenstates are real 
for real-symmetric Hamiltonians. In other words, the state of 
the A^-dimensional dynamical system 

= -diag{Ai,...,AAr}, i(0) = 

must hit a point whose coordinates are integer multiples of tt, 
with the correct parity. As the pariety is not affected by the 
modulo 2tt operation the problem reduces to whether the state 
of the system 

= -diag{Ai,...,AAr} mod 27r, a;(0) = (8) 

hits a point x* with coordinates or tt, depending on the signs 
of the various VikVjk- This dynamical system is the linear 



flow on the iV-torus, and by |j2] Proposition 1.5.1] it is 
minimal |2, Definition 1.3.2], that is, the orbit of every point 
is dense. (Observe that minimality is stronger than topological 
transitivity ||2[ Definition 1.3.1]!) Hence we can get arbitrarily 
close to the point with desired coordinates provided we allow 
t to be large enough. ■ 
To get a quantitative estimate of how close the state x{t) 
has to be to the target point x* with coordinates 0, tt so that 
P{hj) - P{hh t) < e note that 



k 



{±l)VikVjk 



'^VikVjke + VikVjk (±1 

k 



k 

E 

k 



±1 



I 



Thus we have y/p{i,j) ~ \/p{i,j,t) < J2^. |±1 
From physical considerations we know that < p{i,j, t) < 1. 

-P = (VPmax - VP) (VP^ 



Pn 



thus shows that to secure p„ 



p < e, it suffices to make 
|±1 — e^*^*-'^! < 27v- If we set Xk{t) = —Xkt and denote 
the dynamical target state as x^, := or tt, it suffices that the 
target state and the actual state are within the specification 
\xl — Xk(t)\ < sin^^ (2^)- "^ince the topology induced by 
hypercubes is equivalent to the usual topology induced by 
balls, the latter specification can be achieved by the density of 
the orbit of 0. 

B. Estimate of time to attain maximum probability 

The preceding material only tells us that one can reach 
arbitrarily closely the maximum probability, but it does not 
tell us how much time it takes. A conservative estimate can 
be derived, conservative in the sense that it assumes that 
the dynamical evolution in x has been discretized as the 
translation on the torus ||2] Sec. 1.4], 

a;(n + 1) = — diag{Ai, Atv} mod 27r, x(0) = 

The key result is that, under the condition that 1 
are rationally independent, the translation on the torus is also 
minimal |2, Prop. 1.4.1]. In this case, the problem consists in 
finding n e N such that 

-Afcn = (2m + \)tt if Sgiv{vikVjk) = -1 
-Xkn = (2m)7r if Sgn{vikV^j) = +1 

is satisfied with arbitrary accuracy, which can be derived from: 



Theorem 2 (Nowak [14']): For any - 
1, A^, there exist infinitely many {{pi, 
such that 



G R\Z, k = 



N 
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k=l 



A 



Pk 
q 



< 



1 + 1/JV 



(9) 



where the supremum, cjy, of all c^r's satisfying the above 
is known as ci — a/5, C2 = •\/23/2, and for larger N 
estimated as C3 > 1.7739, and, for N > A, cn > {N + 
l)(^+i)/2iV-^/2(^/2)(^+i)/Vr((iV + 5)/2). 

We already know by the minimality of the discrete 
flow, except for the error bound. Therefore, the set of 
((pi, ...jPat), (7) e X N of the above theorem and the 
set of those that satisfy the parity condition have a nonempty 
intersection. Take a {{pi, ...,ppf),q) in this intersection; thus 
the Pk are consistent with the parity condition. Hence 



N 

E 

k=l 



>^kn Pk 



< 



::n 

1+1/N 



Taking n = q yields an £^ error bound of c^^^^/g^/^ on the 



tM 



(j:/7r)-dynamics and furthermore an error bound of 
e on the x-dynamics. Thus the time it takes to be within an 
ball of "radius" e around one of the desired points is estimated 

as n = g « —pf^- 
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